Explicit determination of mean first-passage time for random walks on deterministic 

uniform recursive trees 
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The determination of mean first-passage time (MFPT) for random walks in networks is a theo- 
retical challenge, and is a topic of considerable recent interest within the physics community. In 
this paper, according to the known connections between MFPT, effective resistance, and the eigen- 
values of graph Laplacian, we first study analytically the MFPT between all node pairs of a class of 
growing treelike networks, which we term deterministic uniform recursive trees (DURTs), since one 
of its particular cases is a deterministic version of the famous uniform recursive tree. The interest- 
ing quantity is determined exactly through the recursive relation of the Laplacian spectra obtained 
from the special construction of DURTs. The analytical result shows that the MFPT between all 
couples of nodes in DURTs varies as NlnN for large networks with node number N. Second, we 
study trapping on a particular network of DURTs, focusing on a special case with the immobile 
trap positioned at a node having largest degree. We determine exactly the average trapping time 
(ATT) that is defined as the average of FPT from all nodes to the trap. In contrast to the scaling of 
the MFPT, the leading behavior of ATT is a linear function of N. Interestingly, we show that the 
behavior for ATT of the trapping problem is related to the trapping location, which is in comparison 
with the phenomenon of trapping on fractal T-graph although both networks exhibit tree structure. 
Finally, we believe that the methods could open the way to exactly calculate the MFPT and ATT 
in a wide range of deterministic media. 

PACS numbers: 89.75.Hc, 05.40.Fb, 05.60.Cd, 02.10.Yn 



I. INTRODUCTION 

The field of complex networks has been very active in 
the past decade, since they have been proven a power- 
ful tool to describe very diverse systems in nature and 
society One of the ultimate goals of the study 

of complex networks is to understand the influences of 
network structure on dynamics running on them 
Among various dynamical processes, random walks on 
networks are fundamental to many branches of science 
and engineerin g an d have received a surge of interest in 
recent years [6H13|. As a basic dynamical process, ran- 
dom walks arc relevant to a variety of aspects of complex 
networks, such as target problem [141 . community de- 
tection [ill fl6ll . n etwork routing |l7j. reaction-diffusion 
processesjlcj ll ojj . and so on. Therefore, it is of major 
theoretical interest and practical importance to investi- 
gate random walks on complex networks. 

A primary quantity of interest that relates to random 
walks is the first-passage time (FPT) defined as the ex- 
pected time for a walker to first reach a given destination 
node starting from a source point. The importance of 
FPTs originates from the following main aspects. First, 
their first encounter properties are relevant to those in 
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a plethora of real situations |ll|, such as transport, dis- 
ease spreading, target search, to name a few. On the 
other hand, they can measure the efficiency of random 
navigation on networks I'd . Last but not least, many 
other quantities for random walks can be expressed in 
terms of FPTs, and much information about the random- 
walks dynamics is encoded in FPTs [2l| . Recently, there 
have been a growing number of theoretical studies on 
FPTs H H [ll El. Many authors have devoted their 
endeavors to study the average of FPTs to hit a given 
target node from all other nodes [22M29I ] . In addition, 
relevant work also addressed the FPTs between all cou- 
ples of nodes, giving some numerical results [3(1 Hl| or 
approximate scalings (SrHII, [32-34]. These simulation re- 
sults and scaling laws are necessary as a first step toward 
understanding random walks on networks; however they 
do not provide a complete picture of the random-walk 
dynamics, and analytical exact solutions are helpful in 
this regard [351 ] . 

In this paper, we study analytically random walks on 
a class of deterministic treelike networks. By using the 
links between the random walks, electrical networks, and 
Laplacian spectra, we first compute exactly the mean 
first-passage time (MFPT) between two nodes over all 
pairs of nodes. The obtained explicit formula indicates 
that for large networks with N nodes, the MFPT is 
asymptotic to N In N. Then, we study the trapping prob- 
lem, a particular random- walk issue, on a special case of 
the network family with a trap fixed at a node of the 
highest degree. We derive rigorously the average tapping 
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est [4J-|48| and its relevance to real-life networks. For 
instance, it is small- world [MB IS IIS 51] 5 particularly, 
the so-called border tree motifs have been shown to be 
present, in a significant way, in real- world systems [50| . 
In the rest of this paper, we will study random walks 
performed on DURTs with an aim to better understand 
dynamical process occurring on them. 



III. FORMULATING STANDARD 
RANDOM- WALKS ON DURTS 



FIG. 1: (Color online) The first several iterative processes of 
a particular network for the case of m = 2. 



time (ATT), which is the average of FPTs from all nodes 
to the trap. We show that in contrast to the scaling of 
MFPT, the leading behavior of ATT grows lineally with 
N. Since the MFPT can be considered as the average 
of ATT with the trap distributed uniformly on all nodes 
of the entire network, we conclude that the trap location 
has an important influence on the behavior of ATT. We 
expect that the our analysis technique could be applica- 
ble to determining MFPT and ATT for a broad range of 
deterministic networks. 



II. THE DETERMINISTIC UNIFORM 
RECURSIVE TREES 

We first introduce the model concerned, which is 
a class of trees (networks) defined in an iterative 
method [M]- Let U g (g > 0) denote the networks af- 
ter g iterations. Then the networks can be generated 
as follows. Initially (g = 0), Uq has two nodes con- 
nected by an edge. For g > 1, we can obtain U g from 
Ug—i by adding to (m is a natural number) new off- 
spring nodes to each existing node in U g -\. At each 
iteration gi (gi > I), the number of newly generated 
nodes is L(gi) — 2 m (to + l) 9i_1 . Thus, the order (i.e., 
number of nodes) and the number of edges in U g arc 
N g = Eg,=o L (9i) = 2 (m + iy and K g = N g - 1 = 
2 (to + I) 9 — I, respectively. Figure [T] shows an example 
of the network family for a special case of to = 2 after 
four iterations. 

Note that the network for the special case of m = I 
is actually a deterministic version of the uniform re- 
cursive tree (URT) [371 ] . which is a principal famous 
model [M [M] for random graphs [13] and has a vari- 
ety of important applications in many aspects [4ll - |43j |. 
Moreover, since this particular case of the networks under 
consideration has similar topological characteristics as 
the URT, we call the investigated networks deterministic 
uniform recursive trees (DURTs), which could be helpful 
for better understanding of the nature of the URT. 

We study this model because of its intrinsic inter- 



We study a simple model for random walks on the 
DURTs Ug. At each time step, the walker makes a jump 
from its current location to any of its nearest neighbors 
with uniform probability. We are interested in the FPT 
of a random walker starting from a source to a given 
target point, averaged over all node pairs of source and 
target points. 

To determine the FPT between a pair of two different 
nodes, one can make use of the method of the pseudoin- 
verse of the Laplacian matrix 51| for U g , where random 
walks are performed. It is a very efficient method, which 
allows to obtain the FPT between two arbitrarily dis- 
tinct nodes directly from the network topology and only 
requires inversion of a single N, 
doinversc of the Laplacian matrix 
variant of the inverse of its Laplacian matrix, L s 
elements I 9 , - of the latter are defined as follows: l? 4 
if nodes i and j are connected by a link, otherwise If, 

while 4 ff) = h ( viz., degree of node i). Then, the pseu 
doinverse of the Laplacian matrix L„ is defined to be [52| 



v g x N g matrix. The pseu- 
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where e g is the A^-dimensional "one" vector, i.e., e g = 
(1,1,...,1) T . 

We use Fij(g) to denote the FPT for the walker in U g , 
starting from node i to node j, which can be expressed 
in terms of the entries ijj 9 of as follows [31( : 

F ij (9)=J2(lln 9 -llf-l}': + l}f)l 3 nn , (2) 
n=l 

where Z 9 „ is the n entry of the diagonal of the Laplacian 
matrix h g . So the total, F to t(g), for FPTs between all 
pairs of nodes in U g reads 



(3) 



and the MFPT averaged over all node pairs, (F) g , is then 



(F) B = 
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FIG. 2: (Color online) Mean first-passage time {F) g as a func- 
tion of the iteration g on a semilogarithmic scale for different 
values of parameter m. The filled symbols are the numerical 
results obtained by direct calculation from Eqs. @ and Q, 
while the empty symbols correspond to the analytical values 
from Eq. (|25p . both of which are consistent with each other. 



Equations ([2]) and ((4]) show that the issue of com- 
puting (F) g is reduced to finding the elements of the 
pseudoinverse matrix L g . Since for large g the network 
order N g increases exponentially with g, it becomes in- 
tractable to obtain {F) g through direct calculation using 
the pseudoinverse matrix, because of the limitations of 
time and computer memory, and one can compute di- 
rectly the MFPT only for the first iterations (see Fig. [2j . 
Thus, it would be satisfactory if good methods could be 
proposed to get around this problem. Fortunately, the 
particular construction of the DURTs and the connec- 
tion [53j, |54j between effective resistance and the FPTs 
for random walks allow us to calculate analytically the 
MFPT to obtain a rigorous solution. Details will be pro- 
vided below. 



IV. EXACT SOLUTION TO MFPT AVERAGED 
OVER ALL NODE PAIRS 

In order to avoid the computational complexity of 
inverting the matrix, in what follows, we will use the 
connection between the electronic networks and random 
walks to find a closed-form formula for MFPT, (F) g . 



A. Relation for commute time and effective 
resistance between two nodes 

For a given graph G, its underlying electrical net- 
work [55[ can be obtained by replacing each edge of G 
with a unit resistor. The effective resistance of electrical 
network provides an alternative way to com put e FPTs 
for random walks on the original network [53l. l54j . It has 
been proven that for a connected graph, the FPTs, F^ 
and Fji between nodes i and j, and the effective resis- 
tance, Rij , between these two nodes satisfy the following 



F 



Fji — 2 K Rij , 



(5) 



where K is the number of all edges in the graph and 
Fij is the expected time that a random walker spends on 
reaching node j for the first time, starting from node i. 
Actually, the sum, F^ + Fji, is the average time for a 
walker to go from i to j and back or vice versa, and it is 
often called commute time [56| represented by CV,-, i.e., 
Cij = F^ 



Fji . By symmetry, Cij — Cji . Then, Eq. ([5]) 



can be recast as 



2KRi 



(6) 



Thus, if we view U g as resistor networks by consider- 
ing each edge to be unit resistor, according to the close 
relation between FPTs and effective resistance shown in 
Eq. ([5]) and ©, Eq. © can be rewritten as 



Ftot(g) 



1 N g N g 



(7) 



i^J 3=1 



where Cij(g) and Rij(g) represent respectively the com- 
mute time and effective resistance between two nodes i 
and j of U g . Analogously, Eq. ((4]) can be recast in terms 
of effective resistances as 



Ftot(g) 

N g (N g - 1) 



1 



(8) 



where the sum of effective resistors between all pairs of 
nodes of U g is the so-called Kirchhoff index [57] . which 
we denote by Rtot(g)- Using the previously obtained re- 
sults [H, , the following relation holds: 



Rtot (9)=EE Rii (9) = 2N g J2j- 

.7=1 »=2 1 



(9) 



where (i = 2, . . . , N g ) are all the nonzero eigenvalues 
of Laplacian matrix, L ff , of network U g . Note that since 
Ug is connected, its Laplacian matrix has only one zero 
eigenvalue Ai, i.e., Ai = 0. Then, we have 



N ' 1 



(10) 



Having (F) g in terms of the sum of the reciprocal of all 
nonzero Laplacian eigenvalues, the next step is to deter- 
mine this sum. 



B. Determining MFPT using Laplacian eigenvalues 

After reducing the problem to finding the total of the 
reciprocal of all nonzero eigenvalues of L„, in the follow- 
ing text, we will resolve this problem. 
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By construction, it is easy to derive the following re- 
cursion relation between L g and L g -i, 



order of the Laplacian eigenvalues of U g |48( . Solving Eq. 
(fT6"|) . its two roots are obtained to be 



f L g _i + ml -I -I • • • -I \ 
-I 10 
-I 10 



(11) 



2 ,Arl 



\ -I I / 

where each block is a 2(m + l)*" 1 x 2(m + l) t_1 matrix 
and I is the identity matrix. Then, u sing the elementary 
matrix operations and the results in |6Cf . the character- 
istic polynomial, P g (x) , of L g satisfies [4a | 



and 



v /(Ar 1 +m + l)2-4Ar 1 ), 

(17) 



\ 9 

/v i+2(m+l)s- 



Ar 1 +m+l+ 



Pg(x) = det(xI-Lg) 

= {det^-l)!]}' 



det 



which can be rewritten recursively as 



where 



tp(x) 



(13) 



(14) 



Since there are 2(m + l) 9 nodes in U g , the Laplacian 
matrix L g has 2(m + l) 9 eigenvalues, which are repre- 
sented as Af, A|, . . . , ^2(m+i)s> respectively. We denote 
by Eg the set of these Laplacian eigenvalues, i.e., E g = 
{Af, \ 9 2 , . . . , ^2(m+i)s }' an< ^ we assum e that Af < A?; < 



- A 2(m+1)9- 



The set E g can be classified into two 



subsets represented by Eg and Eg"', respectively [48| . 
That is to say, E g = Eg 1 " 1 U Eg , where Eg 1 " 1 consists of 
eigenvalue 1 with multiplicity 2(m — l)(m + l) 9-1 , 

{1,1,1,..., 1,1}, (15) 



(2) 



2(m-l)(m+l)s- 1 



+ m + l) 2 - 4Af 

(18) 

respectively. Inserting each of the 2(to+1) 9 1 elements of 

into Eqs. (fTTj) and (fT8|) . one obtains the subset _E g 
with cardinality 4(m + 1) 9_1 . Considering Eq — {0,2} 
and recursively applying Eqs. (fTT|) and ([18]). all Laplacian 
eigenvalues of /7 g can be fully determined. 

-vl)%aving obtaining the recursive solutions of the Lapla- 
cian spectra of U g , we continue to calculate the sum on 
the right-hand side of Eq. (fTOf , which is represented by 
S g henceforth. Note that although we fail to determine 
all the eigenvalues of L g in an explicit way, we will show 
that it is possible to provide a closed-form expression for 
S g . By definition, we have 



S 



2(m+l) s ^ 

2J T? 



i=1 



4(m+l)f- 1 3 



(i) » 



4=2 



(19) 



we denote the two sums by Sg 1 ^, and S g 2 \ respectively. 
From Eq. (fT5|) . we can easily get the first sum, 



S^ =2(m-l)(m+l) 9 " 1 . 



(20) 



in which the distinctness of elements has been ignored. 
The remaining 4(m + 1) 9_1 eigenvalues of L g , 

(2) 

forming the subset Eg , are determined by equa- 
tion P g -i((p(x)) = and expressed separately by 
For the sake of convenience, we 



\ 9 \ 9 \ 9 

A x , A 2 , . . . , A 4(m+1)g _ 1 . 

presume Af < Af; < . . . 

\A X , A 2 , . . . , ^4( m+ i)a-i /■ 

According to Eq. (fM] 
bitrary element in E g -\ 
tions of x — m = 



< \ 9 



Thus, E, 



(2) 



it is obvious that for an ar- 
say A 9-1 e E g -i, both solu- 
9-1 belong to E {2) . To fa- 



A 



cilitate the following computation, we rewrite equation 
x — m r -^j = A 9 in an alternative way as 



(Af 1 + m + l)x+ A 



0. 



Moreover, we use notations A? and A 



i+2(m+l)9 



(16) 
_i to repre- 



Thc second sum can be evaluated as 



4(m+l) 



y ~ 




2(m+l)"- 1 7g , Ts 

E A i ~T A i+2(m+l)9-i 
\9\9 



\ y \ y \9 

i=2 * A i+2(m+l)9- 1 A l+2(m+l)0- 1 



(21) 



Since Af and Xf +2 r m+1 \ g -i are the two roots of Eq. (fj"6l) . 
according to the Vieta's formulas, we have Af + 
-i =Af _1 +m+landAfAf +2(m+1)e _ 1 = Af -1 . 



\9 

j+2(m+l)s 



sent it two solutions, which provide a natural increasing Moreover, considering Af = 0, so Af +2(m+1)g _ 1 = m + 1 
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Then, Eq. (|21j) is rewritten as 



2(m+l) fl 



4 2) = £ 



A?" 1 + m + 1 



1 



i=2 



A; 



9-1 



m + 1 

2(m+l)«- 1 



2(m + I) 9 - 1 - 1 + (m + 1) 



1 



1 



A? -1 m + 1 



1 



2(m + l) s - 1 -l + (m + l)S 9 _i + ^ ; . ( 



(22) 



Using J = 5 ff - 4 1} = S g - 2(m - l)(m + l) 9 " 1 and 
after some simplification, Eq. (|22|) becomes 



(23) 



S g = (m+ l)5 9 _i + 2m(m + l) 9 " 1 - 



771 + 1 



With the initial condition So = \ i Eq- (|23| can be solved 
to yield 



(777 + ly 



+ (m + l) 9 - 1 (2mg - 1) + 



Since =2S g ,we have 

= (777 + 1) 9 + 2(777 + 1) 9 " 1 (2777 3 -1) + 



777 +1 



(24) 



(25) 



We have confirmed this closed- form expression for (F) g 
against direct computation from Eqs. ([5J and (j4|). For 
all range of g and different values of m, they completely 
agree with each other, which shows that the analytical 
formula provided by Eq. (|25|) is right. Figure [2] shows the 
comparison between the numerical and predicted results, 
with the latter plotted by the full expression for the sum 
in Eq. {25]). 

We show next how to represent MFPT, (F) g , as a 
function of the network order N g , with the aim to ob- 
tain the relation between these two quantities. Recall- 



ing 



l) 9 , we have (m + l) 9 = and 

9 = 10 gm+i { -t) = lo g m +i N g - log ro+1 2. These rela- 
tions enable one to write (F) g in the following form: 



Ng 
loe 



2(m 
(Ei.) 

-1 \ 2 > 



N„ 



2m 



-N„loe 



777 +1 

2mlog m+1 2 + l 
m + 1 



iV,, 



1 



(26) 



Equation (I26|) unveils the explicit dependence relation 
of MFPT on the network order and parameter 777. For 
large network, i.e., N g — > 00, we have following expres- 
sion: 

(F) 6 



2m 



(777 + 1) ln(?77 + 1) 



N n \nN n 



(27) 



This leading asymptotic N g In N g dependence of MFPT 
with the network order is in contrast to the linear scaling 
previously obtained by numerical simulations for scale- 
free networks, such as the Apollonian networks [30| and 




100 |^ 1000 



FIG. 3: (Color online) Mean first-passage time {F) g versus 
the network order N g on a log-log scale. The filled symbols de- 
scribe the analytical results shown in Eq. (|26[) . The solid lines 
represent the corresponding leading scaling given by Eq. (|27|) . 



the pseudofractal scale- free web j32[. Figure shows how 
the MFPT scales with the network order for two values 
of parameter 777. From Fig. [31 it is clear that for properly 
large network order N g , the dominating term provided 
by Eq. (|2"7| and described by the curve lines agrees well 
with the exact formula given by Eq. (|26)) . 



V. MFPT FOR TRAPPING IN A SPECIAL 
NETWORK 

In the preceding section, we have shown that the 
MFPT averaged over all node pairs, (F) g , varies with 



NghNg. 



Below we 



the network order N g as (F) c 
will show that the scaling for MFPT averaged over part 
of node couples may be different. For this purpose, we 
will study the trapping issue in a particular network for 
777 = 1 case, which is a random- walk problem where a 
trap is positioned at a given location. We focus on a spe- 
cial case with the trap fixed at a node with the largest 
degree (hereafter called hub node) absorbing all particles 
visiting it, which is a simplistic version of trapping in 
complex networks [6lj . 

For simplicity, we continue to use the notation U g to 
represent the network for m — 1 case after g iterations. 
Figure 2] illustrates the first several iterations for this net- 
work. In fact, the network has a self-similar structure, 
which is obvious from the following equivalent construc- 
tion method of the network: suppose one has U g , the next 
generation of the network, f/ 9 +i, can be obtained by join- 
ing two Ug, see Fig. [S] We call the two components, U g , 
in U g +i the original U g and duplicate U g , respectively. 
For the convenience of description, we label all node in 
Ug+x using the following way: the nodes in the original 
U g are labeled as 1 Q , 2 Q , . . ., N 9o , while nodes in the copy 
of U g are labeled as l c , 2 C , . . ., N gc . The trap is located 
at the node belonging to the original U g with label 1 Q . 



6 




g=i g=3 



FIG. 4: (Color online) The growth process for the first three 
generations of the particular network corresponding torn = 1 
case. 



1o 1c 

U, m II — II Uf 



FIG. 5: (Color online) Second construction method of the 
network for m = 1 case, which highlights self-similarity. The 
network after g + 1 generations, U g +i, consists of two replicas 
of U g denoted separately by C/g 1 ' and Ug 2 \ which are con- 
nected to each other by adding a link between two hub nodes 
with one in Ug and the other in Ug 2 \ 



Let Ti{g) be the trapping time (TT) of a node i mU g , 
which is the expected time for a walker starting from % 
to first visit the trap node 1 Q . Obviously, for all g > 0, 
T\ a {g) — 0. We first calculate the quantity ?\ c (<?) that 
is useful for deriving the main result. Since node l c is 
a neighbor of the trap node l c , according to the result 
obtained previously in [62|, i.e., Eq. (9) in (62[, we have 

T lc ( 3 ) = iV g -l = 2 9+1 -l. (28) 

Let Tt t(g) denote the sum of trapping time for all 
nodes in U g , i.e., 

T tot (g)=J2 T i(9)- (29) 

ieu B 

Then, the mean trapping time (MTT), (T) g , which is 
the average of Ti(g) over all initial nodes distributed uni- 
formly in U g , is given by 




Thus, to obtain (T) g , we should first explicitly determine 
the quantity T to t{g), which can be settled using a recur- 
sive way. 

According the second construction method of the net- 
work, it is not difficult to express T to t(g + 1) in terms of 



Ttot(ff)- By definition, we have 

T to t(ff + l) = ]T T 4 (.g + 1)+ T *(.9 + l) 

iet/f 5 *Gt/< 2) 

= £ T *fo)+ ]T [^i c (.9 + l) + F lclo (. 9 + l)] 

ie c/< 2) 

= 2T tot (g) + N g T lc (g + l). (31) 

Considering N g = 29+ 1 , T lc (g + 1) = 2^+ 2 - 1, and the 
initial condition T tot (0) = 1, Eq. (|3T|) is inductively to 
obtain 

Ttot(g) =4 9+1 -2<>(g + 3). (32) 

Plugging the last expression into Eq. (f3T)|). we arrive at 
the closed-form expression for the MTT on network U g 
for the limiting case of m = 1 , 

(T) g = 2» +1 -i±^. (33) 
Clearly, for large network (i.e., N g — > oo), 

(T) g « N g , (34) 

implying that the MTT (T) g increases linearly with the 
network order. 

We have checked the above analytical result using ex- 
tensive simulations. In Fig. [6j we plot the simulation 
results against Eq. (|32|) for different values of g. For all 
values of g, the numerical results are in complete agree- 
ment with the analytical results. Note that the linear de- 
pendence of MTT on network order provided by Eq. (|34l) 
is consistent with the previously obtained results in |63| 
by using a simple approximate method, where it was 
shown, that for the trapping problem in scale-free net- 
works having a degree distribution P(k) ~ fc -7 , when 
the fixed trap is positioned at a node with highest de- 
gree, the MTT (T) varies with the network order N as 
(T) ~ N p with P = (7 - 2)/( 7 - 1). Since for an expo- 
nential network, such as the one addressed here, it can 
be considered as a scale- free network with 7 = 00 [J 0] , 
which leads to /3 = 1, in agreement with the result given 
in Eq. (|34]l . Thus, the exact linear scaling obtained here 
confirms the general case, which was derived based on a 
simple continuous approximation [63j. 

From the above results, we know that the leading be- 
haviors for (F) g and (T) g are evidently different. The 
former follows (F) g ~ N g lnN g , while the latter obeys 
(F) g ~ Ng, less than that of the former. The distinct- 
ness between the two scalings can be interpreted by the 
following heuristic arguments on the basis of the pecu- 
liar structure of the network. In the trapping problem 
addressed here, the location for trap node is particularly 
selected, which lies at a node with largest degree. In fact, 
the trap node is the center of the network (see Fig. [4} : 
one-half of nodes (i.e., descendants of the trap node) lies 
at one side of it, one-half (i.e., descendants of node l c 
including l c itself) at the other side. In this case, the 
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FIG. 6: Mean trapping time {T) g as a function of the itera- 
tion g on a semilogarithmic scale. The filled symbols are the 
numerical results obtained by direct calculation via Eq. 
while the empty symbols correspond to the analytical values 
provided by Eq. (|32[l . The analytical and numerical values 
are consistent with each other. 



walker, irrespective of its starting point, will visit at most 
half region of the whole network before being trapped. 
On the contrary, for some pairs of nodes, such as those 
couples of nodes with both ends being the descendants 
of 1 and l c , respectively, the walker may visit large part 
(even the entire part) of the network before hitting the 
target node. This is the main reason why (T) g is less 
than (F) g . 

Notice that, the random walks discussed in preceding 
section may be considered as a trapping problem with the 
trap uniformly distributed throughout all nodes on the 
networks. The different scalings between (F) g and (T) g 
can lead us to conclude that the location of trap has a 
significant effect on the leading behavior of the MTT for 
trapping problem on one particular network of DURTs, 
which is in sharp comparison with that of trapping de- 
fined on the T graph notwithstanding its tree structure, 
where the MTT is independent of the trap position [rJij . 
The root of this disparity of the behaviors for random 
walks on a DURT and the T graph might lie in their dis- 
tinct structural properties. Although they are both trees, 
the former is small-world with the average distance be- 
having logarithmically with its order [36J ; while the latter 
is not small-world having an average distance increasing 
algebraically with network order [64j |. Particulary, the 
T graph is a fractal, while the DURT is not (its fractal 
dimension is infinite (65j). This fractality has also been 
shown to distinguish diffusion in scale-free networks. For 
details, please see Refs. 0, [28j]. It should be stressed 



that here we only give a possible reason for this differ- 
ence, the genuine explanations need further investigation 
in the future. 

VI. CONCLUSIONS 

We have studied the standard random walks on a fam- 
ily of deterministic treelike networks, exhibiting small- 
world behavior. By applying the connection between the 
FPTs and the Laplacian eigenvalues, we have determined 
explicitly the MFPT averaged over all pairs of nodes in 
the networks. The obtained solution shows that for large 
networks with order N, the MFPT grows approxima- 
tively with N as N\nN. We also presented that com- 
pared to the linear scaling of looped networks, such as the 
Apollonian networks [3(| and the pseudofractal scale-free 
web [13, the DURTs studied here induce a slowing down 
of diffusion dynamics, providing a useful insight into ran- 
dom walks on treelike small- world networks. 

In the second part of this work, we have investigated 
the trapping issue on the a particular network of the 
DURTs, concentrating on a special case with the trap- 
ping positioned at a node with highest degree. We have 
obtained the explicit solution of the ATT, whose lead- 
ing behavior varies lineally with network order. Based 
on the fact that the standard random walks addressed in 
the first part of the work may be looked upon as a general 
trapping problem with trap being distributed uniformly 
on every node in the whole network, we have drawn a 
conclusion that the scaling of ATT for trapping depends 
on the location of trap. Finally, it is expected that the 
analytical computation methods for MFPT and ATT can 
be extended to other deterministic media. 
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